Does the interaction potential determine both the fragility of a liquid and the 
vibrational properties of its glassy state ? 
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By performing molecular dynamics simulations of binary Lennard- Jones systems with three dif- 
ferent potentials, we show that increase of anharmonicity and capacity for intermolecular coupling 
of the potential is the cause of (i) the increase of kinetic fragility and nonexponentiality in the liquid 
state, (ii) the T g -scaled temperature dependence of the nonergodicity parameter determined by the 
vibrations at low temperatures in the glassy state. Naturally these parameters correlate with each 
other, as observed experimentally by T. Scopigno et ai, Science 302, 849 (2003). 
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The structural relaxation time, r, of all glass-forming 
liquids increases on cooling. It becomes so long at some 
temperature T g that equilibrium cannot be maintained 
and the liquid is transformed to a glass. T g is defined 
as the temperature at which r reaches some arbitrarily 
chosen long time, say 10 2 s. Although this behavior is 
shared by glass-formers of diverse chemical and physi- 
cal structures, the scaled temperature dependence of r 
in the liquid state can differ greatly from one liquid to 
another in the degree of departure from the Arrhcnius 
scaled temperature dependence PJQ- The departure can 
be characterized by the rapidity of the change of log(r) 
with Tg/T at T g /T~\, which is given by the steepness 
index or the fragility to defined by Q 
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The values of to of glass-formers of all kinds vary over a 
large range, from the least value of about 17 for strong 
glass-formers (like silica) having nearly Arrhenius scaled 
temperature dependence of t, to values as high as about 
200 found for some glass-formers called fragile. Natu- 
rally, such large variations observed in to beg the ques- 
tion of its microscopic origin. Several attempts have been 
made in the past to correlate to with other dynamic or 
thermodynamic properties, with the hope that the cor- 
relations will lead to the factor or factors that determine 
to. Examples include: (1) The correlation of to with 
n = (1 — (3) at T = T g , where (3 is the stretched ex- 
ponent in the Kohlrausch function, exp[— (i/r) ], used 
to fit the time dependence of the correlation functions 
such as the intermediate scattering functions. (2) The 
correlation of to or (1 — (3) with the mean square dis- 
placement < u 2 > obtained 4] from quasielastic neu- 
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tron scattering measurement of the Debye- Waller factor 
exp[— < u 2 > Q 2 /3}. Glass-former with larger m or 
(1 — (3) has a larger < u 2 > at the same value of T/T g 
and rise more rapidly as a function of T ' /T gi below T g as 
well as near and across T g in the liquid states £j . (3) The 
correlation between to and the slope of the change of the 
configurational entropy, S c , with T/T g at T g jjj. (4) The 
correlation of to with the statistics of potential energy 
minima of the energy landscape 0,0- (5) The correla- 
tion of to with the temperature dependence of the shear 
modulus of the liquid [g . Perhaps the most intriguing of 
all correlations is (6) between to and the vibrational prop- 
erties of the glass at temperatures well below T g found 
recently by T. Scopigno et al. ||. The nonergodicity 
parameter, / (Q,T), at Q=2 nm _1 at very low tempera- 
tures is determined by vibrations. From inelastic X-ray 
scattering data, its temperature dependence is well de- 
scribed by [1 + a (T/Tg)]^ 1 . T. Scopigno showed that to 
and a are proportional for many glass-formers. Appar- 
ently, this last correlation (6) seems to be related to (2) 
for < u 2 {T/T g ) > from neutron scattering at tempera- 
tures well below T g . 

In any glass-former, it is the interaction potential, 
V(r), that determines ultimately all dynamic, thermo- 
dynamic and vibrational properties at all temperatures, 
below and above T g . Changes in an y o f quantities, to, n, 
< u 2 (T/T g ) >, a, S c , free volume v and the degree of 
dynamic heterogeneity, from one glass-former to another 
originate from the change in V(r). Thus, correlations 
found between these quantities are clues for finding out 
which aspects of V(r) determine them and give rise to the 
correlations between them. One would like to examine 
the interaction potentials in real glass- formers. However, 
in such materials, the different kinds of chemical bonding 
and the different sizes of the basic structural unit make 
the comparisons ambiguous. For this reason we consider 
the binary Lennard- Jones particles with different choices 
of interaction potentials V(r) between the particles, and 
perform molecular dynamics (MD) simulations on them 
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to obtain m, n, a, and < u 2 {T/T g ) >. Correlations 
are found between all these quantities, thus reproducing 
the empirical findings from real glass-formers. Since the 
number of particles as well as their density are the same, 
the changes of these quantities are predominantly due to 
the change in V(r). The latter is well controlled, and 
therefore we identify anharmonicity and the capacity of 
intermolecular coupling of V(r) to be responsible for en- 
hancement of m, n, a, and < u 2 (T/T g ) >, and hence 
their correlations. 
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FIG. 1: Potential V(r) governing the A-A interaction. The 
dashed curve is the (8,5) LJ potential for model I, the solid 
curve is the (12,6) LJ potential for model II and the dotted 
curve is the (12,11) LJ potential for model III. 

MD simulations were performed on binary Lennard- 
Jones (LJ) particles systems with three different interac- 
tion potentials by the MD package DLPOLY (TJ. Tech- 
nical details of the MD simulations are given in |l3 • We 
have performed from 10 5 to 6 x 10 7 time steps, depending 
on temperature. We have investigated 12 temperatures 
in the range of [0.675 - 5], [0.416 - 5] and [0.26 - 2] for 
Model I, II and III respectively. All models are com- 
posed of 1500 uncharged particles (1200 species A and 
300 species B). The generalized (q,p) LJ potentials have 
the form, V(r) = (p (^) q -q(f) P ). The param- 

eters ro and Eq represent the position of the minimum 
of the well and its depth, respectively. The reduced LJ 
units 13j are used. The choice of g=12 and p=6 cor- 
responds to the standard LJ potential used by Kob & 
Andersen (K-A) 0] and others for extensive studies by 
simulation. For the purpose of investigating the change 
of dynamics with controlled change of V(r), we developed 
two other models by changing only the exponents, q and 
p, of the LJ potential for the A-A interactions. They are 
(q—8, p=5) and (q=l2, p—11) and shown together with 
the (12,6) LJ potential in Fig.^ The well depth and the 
position of the minimum of V(r) are unchanged, and we 
have kept the standard (12,6) LJ potentials of the K-A 



model for the A - B and B - B interactions, in order to 
retain as much as possible the remarkable ability of the 
K-A model to form a glass upon cooling. The (12,11) 
LJ potential is more harmonic than the classical (12,6) 
LJ potential, while the (8,5) LJ potential is a flat well 
and exceedingly anharmonic. Whenceforth the (12,11), 
(12,6) and (8,5) potentials are referred to as Model I, II 
and III respectively, reminding us that anharmoncity is 
increasing in this order. 

The structure has been briefly studied by the radial 
distribution functions g(r) of the (12,11) and (8,5) mod- 
els, which have been calculated and found very similar to 
that of (12,6) model in the temperature ranges investi- 
gated. Respectively for models I, II and III, the position 
of the first peak of g(r) is at 1.072, 1.066 and 1.057 and 
its width at half the maximum is 0.129, 0.147 and 0.191. 
The g(r) of the three models in the supercooled state is 
similar to that in the liquid regime, ensuring the struc- 
tures have disorder in the mid- and long-range. 

Dynamics have been investigated by computing the 
self Fs(Q,t) and the total F(Q,t) intermediate scatter- 
ing functions of particles A at Qo — 2ir/ro for the three 
models. At high temperatures, Fs(Qo,t) decays linear 
exponentially to zero with a characteristic time of about 

0. 45 close to crossover time t c w 1 to 2 ps used as a 
fundamental time in the Coupling Model (CM) flli fl6|. 
When temperature is lowered, the dynamics slows down 
dramatically and a two-step process appears. This be- 
havior is well described by the Mode Coupling Theory 
(MCT) [l7|. From the F s {Qo,t) at each of these lower 
temperatures, we determine the nonergodicity parameter 
(height of the plateau), fs(Qo>T), the relaxation time, 
ta, and the stretched exponent, (3, from the fit to the sec- 
ond step decay of F s (Qo,t) by f s (Qo,T)exp[-(t/T A f}. 
Shown in Fig. are Fs(Qa,t) versus t/rA of all three 
models at the reference temperature T re f defined by 
TA(T re f)— 46435.8, a very long time in our simulations. 
The values of T ref are 0.688, 0.431 and 0.263 for Models 

1, II and III respectively. 

Shown in Fig. (inset) are j3 of the three models as a 
function of T ref /T. At any T ref /T, (1 - 0) is least for 
model I and largest for model III. At T = T ref , /3=0.69, 

0. 65 and 0.60 respectively for Models I, II and III. Thus 
(1 — f3) increases with anharmonicity. In Fig. [3 log(r^) is 
plotted against T re f/T and the data of the three models 
show systematic change. It can be seen that the slope, 
fragility index to(ta) = d(T^fr) as (. T ref/T) -* 1, in- 
creases monotonically in the order of Models I, II and III. 
The values of m(r^) determined from this latter relation 
are 15.07, 18.57 and 26.58 for Models I, II and III. Al- 
ternatively, the estimated values of m(T A ) based on Sas- 
try's method are 0.195, 0.241 and 0.405 for Models 

1, II and III, respectively. Hence, m(r^) increases with 
anharmonicity. The fragility index found in the present 
study for model II is in good agreement with the value 
determined in an earlier work on the same model Q. 
The diffusion coefficient, D A , of particles A was calcu- 
lated from the mean-square displacement < u 2 {t) > at 
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FIG. 2: Self intermediate scattering function Fs(Qo, t) versus 
scaled time t/TA- Dashed, continuous and dotted lines are 
for Models I, II and III respectively. For all three models, 
TA(T r ef)= 46435.8. The inset shows the stretched exponent 
/3 — (1 — n) as a function of the scaled reciprocal temperature 
T ref /T for the three models. (♦) Model I, (•) Model II and 
(■) Model III. For definitions of ta and /?, see text. 



long times when < u 2 (t) > assumes the linear t depen- 
dence. Here another T' re * is defined as the temperature 
at which Da is equal to 1.86 x 10 -5 . In the inset of Fig.|3 
we plot \og{\/ Da) against T' re ^jT for each of the three 
models. They exhibit the same pattern as log(-r>i). Again 



__ dlog(l/D A ) 



as 



the steepness or fragility index, m(DA) — -gpp — jm 

(T£ e f/T) — » 1, increases monotonically in the order of 
Models I, II and III, or with anharmonicity. 

So far we are concerned for dynamic quantities and 
their correlations for T > T re f, as analogues of them in 
the liquid state of real glass-formers. Next we examine 
the vibrational properties of the three models at temper- 
atures lower and much lower than T re f. At low temper- 
atures, relaxation of any kind is absent in the simulation 
time window, and the nonergodicity parameter f(Q,T) 
determined from F(Q,t) is contributed entirely from vi- 
brations. As performed in 0, we have followed the be- 
haviour of f(Q — > 0,T) which has been obtained by a 
Q-quadratic extrapolation of f(Q, T) for the lowest tem- 
peratures. The results of the three models are shown by 
a plot of f(Q -> 0, T)- 1 versus T/T ref in Fig. H 

In all three cases, the dependence of /(Q,T) _1 on 
T/T re f is approximately linear, and f(Q,T)~ 1 has the 
extrapolated value of unity at the origin. The depen- 
dence of f(Q,T) from simulation on T/T re f is governed 
by the parameter, a, through the expression, 



f(Q,TY 



= 1 



T 



T, 



ref 
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just like a similar expression used to represent the depen- 
dence of f(Q,T) on T/Tg of real glass- formers obtained 



FIG. 3: The relaxation times ta obtained from Fs(Qo,t) for 
the three models as a function of T re //T where T re f is defined 
as the temperature at which ta reaches 46435.8. (♦) Model 
I, (•) Model II and (■) Model III. T re f is the analogue of T g 
for simulations when the dynamics of the system slows down 
to more than 10 ns. In the inset, the reciprocal of diffusivity 
Da of A particles for the three models are given as a function 
of T' ra cjT. Here T' trS is now defined as the temperature at 
which D A is equal to 1.857432 x 10" 5 . 
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FIG. 4: /(Q,T) _1 versus T/T ref for the three models. (♦) 
Model I, (•) Model II and (■) Model III. f(Q, T) _1 is almost 
linear relative to T/T re f with a slope noted a. The inset 
shows the correlation of the fragility m with a from the results 
of the three models. 



by inelastic X-ray scattering We see from Fig. that 
the increase of f(Q, T)~ x with T/T re f is fastest for model 
III and slowest for model I. Equivalently stated, the slope 
a is largest for model III and smallest for model I. a in- 
creases with anharmonicity of the potential, like m{rA) or 
to(Da), and (1 — 0), as seen before (Figs. ^ 13 EJ • Hence 
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the interacting potential is the origin of the correlation of 
a with iti(ta) or to(LLi), and (1 — 0) in our simulation, 
suggesting the same holds for real glass-formers. More- 
over, in the inset of Fig.0] we observe that m(Tyi) and a 
are proportional together as shown in |9j . 

Presumably there is no disagreement that the interac- 
tion potential is pivotal in determining all dynamic, ther- 
modynamic and vibrational properties of glass-formers at 
all temperatures. With all these experimentally accessi- 
ble properties originating from the interaction potential, 
it is not surprising to find correlations or anticorrclations 
between them. Making this point is one of the motiva- 
tions of the work in demonstrating that the various corre- 
lations between a, to, (1— /3) and < u 2 (T/T g ) > observed 
in real glass-formers are reproduced by simulations as the 
analogues of correlations between a, m (ta), to (Da) and 
(1 — (3) by varying the interaction potential V(r). The 
results confer a bonus in identifying which feature of the 
interaction potential is responsible for enhancement of a, 
to (ta), m (Da) and (1 — (3). Certainly the anharmonic- 
ity of V(r) increases when going from Model I to II and 
III, but one also can observe from Fig. ^ that the en- 
ergy barrier becomes lower and flatter. The latter trend 
means that neighboring LJ particles are more coupled 
in their motions. The increase in interparticle coupling 
from Model I to Model III is consistent with the posi- 
tion of the first peak of radial distribution functions g(r) 
(1.072, 1.066 and 1.057 for Models I, II and III respec- 
tively) and with its width at half the maximum (0.129, 
0.147 and 0.191 for Models I, II and III). This insight 
from simulation, when transferred to real glass- formers, 
suggests that the capacity for intermolecular coupling 
and anharmonicity of the interaction potential determine 
the dynamic, thermodynamic and vibrational properties 
of glass- formers above as well as below T g . Although 
the thermodynamic variables, configurational entropy S c 
and free volume v, are determined by V(r), they reenter 
into the dynamics by their influence on molecular mo- 
bility. Thus two factors govern dynamics, the capacity 
for intermolecular coupling directly from V(r), and S c 
and v that come indirectly through V(r). On super- 
cooling a liquid, S c and v change, and since the kinetic 
fragility m is the slope of T^-scaled temperature varia- 
tion of t, it is unsurprising that m is correlated with the 
slope of the corresponding change of S c , i.e., the thermo- 
dynamic fragility. The capacity for intermolecular cou- 
pling of V(r) is solely responsible for the shape of the 



dispersion or the nonexponentiality parameter (1 — (3), 
and it also determines r in conjunction with S c and v. 
The results of our simulation with the three potentials 
support this view. Increasing the density of the parti- 
cles of the binary Lennard-Jones system with the fixed 
(12,6) potential effectively forces the particles to be closer 
to each other and thereby increases intermolecular cou- 
pling. The simulation performed in this manner 7] show- 
ing that to increases with density can be reinterpreted as 
due to increase of intermolecular coupling. Intermolecu- 
lar coupling manifests itself in the dynamic properties in 
various ways. The most direct way is the width of the 
dispersion measured by (1 — (3). The others include : (1) 
the Q -2 ^ 1- ' 3 ' dependence of r obtained by quasielastic 
neutron scattering £| , and (2) the proportionality be- 
tween (1 — (3) and (log(r) — log (173)) at constant log(r) or 
at T g (where Tp is the Johari-Goldstein relaxation time). 
All these are evidences for intermolecular coupling that 
must be taken into consideration in conjunction with S c 
and v for explaining all observed experimental facts in 
the liquid state. At low temperatures and deep in the 
glassy state, S c and v having constant values cannot in- 
fluence the temperature dependence of the vibrational 
properties characterized by a. Hence, a is controlled by 
the anharmonicity of V(r), as demonstrated by the sim- 
ulations. 

In summary, we demonstrate by using three different 
interparticle potentials of binary Lennard-Jones systems 
that the capacity for intermolecular coupling and anhar- 
monicity of the potential are responsible for the correla- 
tions between various dynamic, thermodynamic and vi- 
brational properties of glass-formers. Increase of the ca- 
pacity for intermolecular coupling and anharmonicity has 
the effects of increasing the kinetic fragility, to, and the 
nonexponentiality parameter, (1 — f3), in the liquid state, 
and of increasing in the glassy state the parameter a 
that characterize the T g -scaled temperature dependence 
of the nonergodicity parameter determined by vibrations 
at low temperatures. The correlations between to, (1— (3), 
a and other quantities follow as consequences, and their 
observations by experiments explained. 
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